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Abstract 

We investigated the effect due to perpendicular magnetic field on quantum wires where spin- 
orbit interaction (SOI) of electrons is not neglected. Based on the calculated energy dispersion, the 
nonlinear ballistic conductance [G] and electron-diffusion thermoelectric power (S^) are calculated 
as functions of electron density, temperature and applied bias voltage. A low-temperature half-step 
feature in G, which was observed experimentally by Quay et al. [see Nature Physics 6, 336 (2010)], 
as well as a new peak in Sd are reproduced here in the low density regime. These phenomena 
are related to the occurrence of the Zeeman splitting and SOI induced saddle point in the band 
structure, where the channel chemical potential lies within an anticrossing gap between the saddle 
point of the lower subband and the bottom of the upper subband. Additionally, side peaks in G far 
away from the zero bias for the nonlinear transport, as well as a quadratic bias- voltage dependence 
of G near zero voltage, are predicted and discussed. 
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I. INTRODUCTION 



Quasi-one- dimensional (quasi-lD) electron systems (also called quantum wires) have at- 
tracted considerable attention both theoretically-"- and experimentally-'^ since their ex- 
perimental realization in 1986.- When an external magnetic field is applied in a direction 
perpendicular to a quantum wire, Landau quantization in such a system combines with 
lateral confinement, enhancing the quantum-wire confinement potential as well as the ef- 
fective mass of electrons along the wire direction at the same time.-*^ Several studies were 
reported for such systems, including magneto-optical properties-i^ and magneto-transport 
properties—"— in the diffusive regime. Moreover, the ballistic transport properties—"— of 
quantum wires were also explored under a tilted magnetic field. 

On the other hand, the effect of Rashba spin-orbit interaction (SOI) on the ballistic 
and diffusive electron transport properties in quantum wires has also been investigated 
theoretically and experimentally in the absence of a magnetic field.—"— Compared with the 
case when SOI is neglected,— the quantum-wire conductance quantization in the presence 
of SOI is modified significantly at low temperatures.— Although the SOI-induced gap has 
been observed experimentally,™ there are only a few theoretical studies dealing with ballistic 
electron transport and thermoelectric properties in parabolically-confined quantum wires 
with both SOI and a magnetic field.—"— 

In this paper, we include the effects due to the edges of a realistic quantum wire through 
sharp and high potential walls forming the boundaries, instead of assuming a parabolic 
potential for the lateral confinement of the quantum- wire . Extending our previous study for 
zero magnetic field in the linear response regime,— we calculate the bias- voltage dependent 
nonlinear ballistic conductance and electron-diffusion thermoelectric power for this quasi- 
one-dimensional system with a channel length less than the electron mean free path, where we 
have included the effect on the SOI of a finite magnetic field perpendicular to the quantum 
wire as well as the Zeeman effect on the electron spin states. A half-step feature in the 
quantized conductance is obtained for low electron densities and low temperatures. 

The rest of this paper is arranged as follows. In SeclTTl we present our model for quantum 
wires with SOI under a perpendicular magnetic field. Numerical results and discussion of 
the energy bands, conductance and diffusion thermoelectric power as functions of electron 
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density, temperature and bias voltage are presented in Sec.lIIII We summarize our results 
in Sec.IVl 



II. THE MODEL 



In the presence of a magnetic field B (in the z direction) perpendicular to a quantum 
wire which is confined in the x direction such that |x| < VV/2 and extends along the 
y direction, the spin-coupled Schrodinger equations for the electron wave-function spinor, 
{x\ ipj. (x)]'^, with the Landau gauge are written as— 
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where /3 and a are the SOI parameters in the transverse x and z directions, respectively. 
Also, e is the eigenenergy of the electrons, the symbols ^ and J, denote the up and down spin 
states, ky is the electron wave number along the wire, m* is the effective mass, Uc = eB/m* 
is the cyclotron frequency, £c ■ 



'h/eB is the magnetic length, and g* is the effective Lande 
spin factor. Additionally, the quantum confinement factor— introduced in the above 
equations for a non-parabolic confinement is given by 
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(3) 



where £o describes the steepness of the confinement potential at the two edges of the wire. 
Compared with the case of 5 = O,-^ we see that ky is replaced by {ky + well as the 

appearance of a Zeeman splitting term, in Eqs. ([1]) and ([2]) when B ^ 

For ballistic electrons in a quantum wire, the charge [n = 0) and the heat (n = 1) currents 
are calculated as— 



3 



^ — oo 

X 'e{-v,,^) ff^l (i - /S) + (i - 
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In this notation, j = 1, 2, ■ ■ ■ is a composite index labehng the spin-resolved subbands, 
Vjky = (1/^) dejky/dky is the group velocity of electrons in the jth subband, 9{x) is the unit 
step function, /Xc is the chemical potential of the quantum wire, /j^^ = {1 + exp[{ejky — 
fis)/kBT]}~^, T is the temperature of the system, s = R, L denotes the right and left 
electrodes connected to two ends of a quantum wire in the y direction, /il = /^c — eVb/2 and 
/^R = A^c + eVb/2 represent the chemical potentials of the left and right electrodes, satisfying 
f^R — f^L = eVfe with Vb being the external bias voltage. For fixed chosen linear electron 
density uid and T, fi^ can be calculated using the following constraint equation 
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i^jky - IJ'c)/kBT 

Once the Q^"\Vb) in Eq. (H]) has been calculated, the bias-dependent ballistic conductance 
G{Vb) of the quantum wire is simply given by 



G{Vb) 



Vb 



(6) 



and the diffusion thermoelectric power S^iVb) may be calculated according to 



Q^'KVb) 



We now employ these quantities to carry out our numerical calculations. 



(7) 



III. NUMERICAL RESULTS 



In our calculations, we define two spin-related characteristic lengths, namely, ia = 
h'^/{2m*a) and £^ = / (2m* f3) . Consequently, four corresponding dimensionless quan- 
tities related to spin-orbit coupling, quantum confinement factor and magnetic-field effect 
are given by: = >V/4 oc a, r^j = W/if^ (x (3, Tq = W/io and n = W/4 oc y/B. The 
scaled electron kinetic energy is Ejky/Eo with Eq = n'^h^ /2m*W'^, the scaled electron wave 
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number is kyW, and the scaled position variable is x/W. Parameters in our calculations 
were chosen as m* = 0.067 mo (as in GaAs) with mo being the free-electron mass, g* = —0.44 
due to exchange interaction, = 10, W = 568. TA and tq = 10^. Other parameters, such 
as r„, B, T, and riio, will be given in the figure captions. 

Figured] shows the calculated eigenenergy Ejky/Eo as a function of kyW for three different 
cases, i.e., (a) r„ = 10, 5 = OT; (b) r„ = 10, 5 = 1 T; (c) r„ = 0, 5 = 1 T. From Fig.f^a), 
we find that the spin-orbit interaction under a ^ couples spin-up and spin-down electron 
states in the absence of a -B field. Here, the energy dispersion curves of the spin-up and 
spin-down states are shifted relative to each other with respect to ky = 0. In the presence 
of a -B field, as shown in Fig.lUJ^b), many anticrossing gaps are created due to both Zeeman 
splitting and SOI, the effective mass of electrons along the wire is enhanced, and new saddle 
points are formed. After the spin-orbit interaction is turned off, we simply see a set of 
Zeeman split parabolas, as demonstrated by Fig.[T]^c). When i? — ?► oo, these subbands will 
convert into discrete Landau levels (with infinite effective mass along the y direction). 

The results for the quantum- wire chemical potential (Ic/Eq, the ballistic conductance 
G [in units of (e'^/h)], and the diffusion thermoelectric power Sd [in units of {—ks/e)] are 
compared in Figs.lU^a), (b) and (c), respectively, for the three cases discussed in Fig.[TJ Here, 
T = 4K and a small bias voltage Vb = 0.1 meV is used for the linear-response regime. From 
Fig.[2]^a) we see that when a 7^ 0, yUc is reduced at high electron densities but enhanced 
at low densities by the B field due to the formed anticrossing gap and increased electron 
effective mass. Additionally, there exists a series of kinks in fic as a function of riiD, which 
is related to the onset of a new subband population. By turning off the SOI, is increased 
for finite values of B and the kinks in fi^ are shifted at the same time due to the suppression 
of the fiattened region around ky = where the density of states is expected to be high. 
This series of kinks observed in Fig.[2]^a) leads to the occurrence of quantum steps for G in 
Fig.[2]^b). Furthermore, each quantum step in G relates to a peak in S^, as presented in 
Fig.[2]^c). It is interesting to note that the application of the B field has a very profound 
effect on the spin-orbit interaction of electrons at low densities, where a new step (indicated 
by a black arrow) in (b) and a new peak in (c) are seen in correspondence to fic lying inside 
the lowest anticrossing gap. Similar half-step features were observed experimentally.— 

Comparisons for /ic, G and Sd as a function of temperature T for the three cases discussed 
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in Fig.[T]are presented in Figs.[3]^a), (b) and (c), respectively, where we chose niD = lO^cm"^ 
and Vb = 0.1 meV. With a finite SOI, the B field reduces (increases) fic in Fig-EI^a) for 
T < 4K (T > 4K). For a = 0, /ic is significantly enhanced. Meanwhile, /ic, which lies 
within the gap of the second and third parabolas, increases with T in contrast to the case 
when a ^ where a decreasing fi^ with T is obtained due to the fact that /ic stays very 
close to the two minima of the third subband in Figs.[T]^a) and (b). For nonzero SOI, G 
increases with T in Fig.|3]^b) as a result of the gradual thermal population of higher subbands 
slightly above /ic. On the other hand, G is decreased when T is increased in the absence of 
spin-orbit couphng since G oc [—dfjl^^/dejky] whose peak strength is reduced by increasing 
T provided /ic lies within a gap. In addition, as can be seen from Fig.lHl^c), monotonically 
increases with T, and becomes linearly proportional to T above 4K.— However, it becomes 
exponentially small for T ^ 1 K. For a = 0, becomes much smaller than in the previous 
two cases due to a larger value of G. 

We present in Fig. |4] comparisons between the bias voltage (Vf,) dependence of G in (a), 
the low- voltage blow-out view of Fig.Hl^a) in (b) and the charge current Q^^-* in (c) for the 
three cases considered in Fig.[Tl where uid = lO^cm"^ and T = 4K. When = 10 in 
Fig.m^a), G initially increases with \Vb\ from 2e^/h in a symmetric way, which is followed 
by a decrease before it increases again. The B field does not alter the shape of the G — Vb 
curve except for a small reduction of G. On the other hand, G decreases with \Vb\ from 
Ae'^/h and eventually saturates at 3e'^/h for Tq = 0. These dramatically distinct behaviors 
in G are attributed to different positions of the chemical potentials in the two electrodes 
relative to the bottom of subbands as V;, — )■ 0, as explained in Fig.O In order to see 
the behavior of G in the low-voltage regime, we display in Fig.Hl^b) the blow-out view of 
Fig. 111(a), where a quadratic dependence of G on H can be identified. Figure Hl^c) exhibits 
a different dependence of Q*^*^-* on Vb- For the cases with Tq, = 10, the leading term in the 
charge current is Q^^^ oc Vb [—dfjky/dejky] for \Vb\ < 1.5 meV. However, for \Vb\ > 2.5 meV in 
Fig.m^c), we find that the leading term changes to Q^^^ oc V^[—d^ fj'^^/de^i^^]. This explains 
the observations in Fig.Hl^b), i.e., G = 2e^/h for VJ, — (linear response) and G oc Vj,^ 
(nonlinear response) for large values of |Vft|. As Tq, = 0, the slope in the linear- response 
regime is greatly enhanced because of a larger value of the Fermi velocity (large /ic), leading 
to a quadratic decrease of G with Vb in Fig.|l](b) as IHI increases. 
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In Fig.|5]we display two-dimensional (2D) contour plots of G(niD, T) in (a) and (b) as 
well as 2D contour plots of S'd(r;,iD, T) in (c) and (d) when B = [(a) and (c)] and when 
-B = 1 T [(b) and (d)]. We chose VJ, = 0.1 meV and = 10. In general, for both B = and 
i? = 1 T, we clearly see quantized steps in G and peaks in with increasing at any 
value of T in this figure. The step and peak positions are almost independent of T. However, 
in the low-density regime when niD < 0.5 x lO^cm"^, the presence of the B field modifies 
the spin-orbit coupling and dramatically alters the behaviors of G and as functions of 
both niD and T. 

We compare the 2D contour plots of G{T, Vb) in Fig.[6]at i? = in (a) and i? = 1 T in 
(b), for chosen niD = lO^cm"^ and = 10. It is easy to see that the two strong peaks 
in G on either side of VJ, = are greatly suppressed by the magnetic field, which is further 
accompanied by an expanded dip around Vf, = 0. In addition, the peaks in G decrease with 
T. 

For completeness, we present in Fig.[7]the 2D contour plots of G(?2id, H) for 5 = and 
i? = 1 T in (a) and (b), respectively, at T = 4K when = 10. It is clear that the quantized 
steps in G become most observable close to H = 0. The B field does not significantly affect 
G in the high-density regime (nio > 2.5 x 10^ cm^^). In the intermediate-density regime 
(1.2 X 10® cm^^ < riiD < 2.5 x 10® cm^^), G only weakly depends on | Vf,| and shows no -B-field 
effect. In the low-density regime (niD < 1.2 x 10® cm"^), however, the B field changes G 
substantially as functions of both niD and \Vb\. The dip around Vb = evolves into a peak 
and vice versa with increasing riiD due to the variation of the relative position of /ic with 
respect to the bottom of the subband. 

IV. SUMMARY AND CONCLUDING REMARKS 

In summary, we have investigated the combined effect of SOI and an external magnetic 
field which is applied perpendicularly to a quantum wire. The quantities calculated are 
the the quasi-one-dimensional energy subbands, the nonlinear ballistic conductance (G) and 
electron-diffusion thermoelectric power (Sd) as functions of electron density, temperature and 
applied bias voltage. An important outcome of our numerical calculations is the observation 
of a half-step in the low-temperature quantized conductance at low electron densities. We 
have accounted for this feature as being due to the occurrence of a saddle point in the 
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subband energy dispersion relation. The saddle point is due to the combined action of the 
applied magnetic field and spin-orbit coupling in the wire. A similar low-temperature half- 
step feature in G was observed experimentally by Quay et al.—. A new peak in Sa is produced 
in the low density regime, and is also related to the occurrence of the magnetic field and 
SOI induced saddle point in the band structure, where the channel chemical potential lies 
within an anticrossing gap between the saddle point of the lower subband and the bottom 
of the upper subband. Additionally, the present nonlinear transport theory has predicted 
side peaks in the quantized conductance as a function of bias voltage away from zero and a 
quadratic bias dependence near zero voltage. This extension of our previous work^'^' clearly 
demonstrates that the application of magnetic field leads to observable differences in the 
ballistic transport properties of a quantum wire structure with SOI playing a role. 
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FIG. 1: (Color online) The calculated eigenenergies ejky/EQ as functions of kyW for the three cases 
considered: (a) Tq, = 10 and = OT; (b) Tq, = 10 and -B = 1 T; (c) Tq, = and B = IT. The 
Zeeman splitting in (c) can be resolved for B = IT. The degeneracy of t and J, spin states at 
ky = are split by the Zeeman effect in (b) to form an anticrossing gap. 
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FIG. 2: (Color online) (a) Chemical potential [i^ (in units of £^0) for a quantum wire, (b) ballistic 
conductance G (in units of {e^ /h)) and (c) diffusion thermoelectric power 5^ (in units of {—kB/e)) 
as functions of the linear density niD for the three = 10, S = T (blue solid curves); 

Tq, = 10, B = IT (red solid curves); and Tq = 0, i? = 1 T (black dashed curves), where T = 4K 
and Vh = 0.1 meV. The black arrow in (b) indicates a half-step in quantized conductance G. 
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FIG. 3: (Color online) Calculated (a) /^cj (b) G and (c) S^, as functions of temperature T for the 
three cases: Tq, = 10, i? = T (blue solid curves); = 10, -B = 1 T (red solid curves); and = 0, 
= 1 T (black dashed curves), where nio = 10^ cm~^ and 14 = 0.1 meV. 
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FIG. 4: (Color online) Calculated (a) G, (b) its blow-out view and (c) charge current Q^*^^ as 
functions of the bias voltage Vf, for = 10, S = OT (blue solid curves), Tq, = 10, i? = 1 T (red 
solid curves), and Tq, = 0, S = 1 T (black dashed curves). Here, T = 4K and tiid = lO^cm"^. 
The vertical dotted line in (b) indicates the position for Vf, = 0. 



13 




FIG. 5: (Color online) Contour plots of G [in (a) and (b)] and Sd [in (c) and (d)] as functions of 
niD (horizontal axis) and T (vertical axis) for B = OT, in ((a) and (c)) and B = IT ((b) and (d)). 
Here, Vf, = 0.1 meV and = 10. The color scale of G in (a) and (b) varies from (blue) to 6.5 
(red), whereas that of Sd in (c) and (d) varies from (blue) to 0.75 (red). 
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FIG. 6: (Color online) Contour plots of G as functions of T (horizontal axis) and Vb (vertical axis) 
for (a) S = T and (b) B = IT. We chose niD = 10^ cm^^ and Tq, = 10. The color scale for G in 
(a) and (b) varies from 1.9 (blue) to 2.9 (red). 




FIG. 7: (Color online) Contour plots of G as functions of nio (transverse) and Vf, (vertical) for 
B = QT [in (a)] and B = IT [in (b)]. Here, T = 4K and Tq = 10. The color scale of G in (a) and 
(b) is from 0.7 (blue) to 6.5 (red). 
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